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ABSTRACT 

Usual treatment of the Nambu- Jona-Lasinio (NJL) model using loop momen- 
tum cutoff suffers from ambiguities in choosing the loop momenta to be cut off 
and violation of (external) gauge invariance. We define the NJL model from the 
starting Lagrangian level by using a higher derivative fermion kinetic term with a 
cutoff parameter A. This definition is free from such ambiguities and manifestly 
keeps the chiral symmetry as well as the gauge invariance. Quantization of this 
higher derivative system, current operators and calculational method are discussed 
in some detail. Calculating the pion decay constant and 7r° — > 27 decay ampli- 
tude, we explicitly demonstrate that the low energy theorem holds. It is observed 
that the NJL mass relation m a = 2mo between the a meson and fermion masses 
no longer holds in this model. We also present a simplified calculational method 
which is valid when the cutoff parameter A is much larger than the energy scale of 
physics. 



1. Introduction 



The Nambu-Jona-Lasinio (NJL) model ' is not denned without an ultraviolet 
cutoff A. Usually this cutoff is introduced at the loop graph level. For instance, 
in the most popular treatment, the fermion one- loop integral is first rewritten 
using the Feynman parameter formula and the ultraviolet cutoff is made for the 
loop momentum variable with which the denominator of the integrand becomes an 
even function. Another treatment, which was adopted in the original NJL paper) 1 ' 
utilizes the dispersion relation. In either treatment, the cutoff is made graph by 
graph (or for each Green function separately). This implies that the theory is not 
defined at the starting Lagrangian level. This is very unsatisfactory. 

Other serious problems in the usual treatment are the ambiguity in introducing 
the cutoff and the consistency with symmetries. The values of the divergent loop 
diagrams depend on the choice of the loop momentum variables for which we make 
the ultraviolet cutoff. We should specify the choice procedure unambiguously. 
Moreover the procedure has to be shown to be consistent with the chiral symmetry 
at least. To show this consistency would not be an easy task if the cutoff is 
introduced graph by graph. Furthermore if the system is coupled to external gauge 
fields the cutoff procedure should also be consistent with the gauge invariance. It 
is indeed very difficult to satisfy the gauge invariance if we adopt such a simple 
cutoff for the loop momentum variables. 

A best example for the last problem is given by vector 2-point functions. Con- 
sider a NJL model which contains 4-fermi interaction in the vector channel and 
couples to the external photon field A^. If one rewrite the 4-fermi interaction term 
by introducing a vector auxiliary field V^, the fermion kinetic term gets to take the 
form 

^{d^-iV^-ieA^. (1.1) 

Note the pararellism between the vector field and the photon field A^. In 
calculating the fermion one-loop diagrams, there is no difference between them; 
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namely, exactly the same Feynman diagrams appear for the vector-vector, vector- 
photon and photon-photon 2-point functions. If we adopt a loop momentum cutoff 
to the diagram, we would obtain gauge non-invariant function not proportional to 
9{ivP 2 —p^P h ' i which may be a good result for the vector-vector function but is clearly 
unacceptable for the photon-photon and vector-photon functions. To achieve the 
gauge invariance for the latter, one sometimes throw away the gauge non-invariant 
piece, quadratic mass term oc g^A 2 , by hand. This might be a correct procedure. 
But if so, then, should we do the same also for the vector-vector case? 

We propose in this paper to define the NJL model by using a higher derivative 
kinetic term for the fermion. We replace the usual fermion kinetic term by the 
following higher derivative one: 

+ • (1.2) 

Here the parameter A 2 plays the role of an ultraviolet cutoff. Note that this effective 
cutoff is made on each fermion propagator but not on each loop momentum, and 
that any loop diagrams are now well-defined and independent of the choice of the 
loop momenta. It should be emphasized that this defines the NJL model already at 
the starting Lagrangian level. It is manifest that this higher derivative definition 
keeps the important chiral symmetry of the system. Moreover if we switch on 
external gauge interactions (such as weak-electromagnetic ones), we can keep the 
gauge-invariance also by simply replacing each derivative factor <f) by a covariant 
one Tp. 

As a mere regularization method there is dimensional regularization or zeta 
function regularization which respects the gauge invariance. But what we need 
here is not a mere regularization but a definition of the NJL model by a 'regular- 
ization' which is not removed forever. If we define the NJL model by dimensional 
regularization, the model will be defined in 4 — e dimensions with a certain e fixed 
and then the physical meaning would become very unclear. In our definition by 
higher derivative, the cutoff parameter A has a physical meaning as the energy 
scale which gives an upper limit of the applicability of the model. 
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This paper is organized as follows. In Sect. 2 we set up the NJL model which we 
discuss in this paper. Since the model is a higher derivative system rather different 
from the canonical one, we discuss there the problems of quantization, current 
operators and coupling to the external gauge fields in detail. In Sect. 3 we present 
a method for reducing the calculations of Feynman diagrams in this system to 
those of the usual first-order derivative fermion system. We calculate some Green 
functions and note, in particular that the scalar meson has a pole below the two 
fermion threshold. In Sect. 4 n° — > 27 amplitude is calculated and is explicitly 
demonstrated to be consistent with the low energy theorem. When the cutoff A is 
much larger than the energy scale we discuss, the calculation can be made much 
simpler. This is shown in Sect. 5. Sect. 6 is devoted to conclusion. In the Appendix, 
we present a generalized Noether procedure of constructing conserved current in a 
generic higher derivative system. 

2. NJL Model Defined by Higher Derivative Kinetic Term 

2.1. NJL Model and QCD-analogue Model 

We consider a fermion tp = (ipif) which carries SU(N C ) color index % and 
SU(ri{) flavor index /. The fundamental representation matrices are denoted by 
T A (A = 1,---,A" C 2 - 1) for the color SU(N C ) and by A a (a = 1, • • • , n\ - 1) 
for the flavor SU(rif), respectively. They are normalized by ti(T A T B ) = \8 AB , 
tr(A a A b ) = 25 ab . We also use flavor singlet matrix A = \j2~jnt l nf proportional to 
rif x rif unit matrix l nt . The NJL model we consider in this paper is the following 
system which possesses chiral U(rif)n x U(n{)i, symmetry up to the axial U(1)a 
anomaly: 

+ IP { [(^ A » 2 + (#T5AV) 2 ] + [(W) 2 + (^75A>) 2 ] } 

" fff^AV) 2 + (VwsAV) 2 ] - IpV^A^) 2 + (^75A>) 2 ] . 

(2.1) 
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Here all fermion bilinears are of color- singlet. The fermion kinetic term is taken to 
be a higher derivative one as explained in the Introduction. This fully defines the 
model from the starting Lagrangian. 

If we introduce auxiliary fields following the well-known procedure, this La- 
grangian can equivalently be rewritten into 



c = i> 



iP{l + W)-M]*-^ (Et E ) + || tr (R»R» + L%J 



+ 



A 2 
A 2 



4 V3vo 9y 



tr(i? M A°)] 2 +[tr(L M A°)] : 



(2.2) 



M(x) = X(x)V R + T$(x)V L ~ -fRp(x)V R - -fLp(x)V L ■ 

with chiral projection operators Vr = (1 + 7s)/2 and Pl = (1 — 7s)/2. Eq.(2.2) 
gives equations of motion for the auxiliary fields as follows: 



£ = 



-R^ — 



■ || [A^rA>l) + A W> L )] 



2A 2 

A 
2A 2 



^A°(^r7mAVr) + A a (<M7„A>R) 
L 9v 

^A°(^l7mAVl) + A a (^L7MA>L) 
9v 



(2.3) 



with V'RjL = PRjL.V'- Note that all the auxiliary fields are rif x rif flavor matrices; 
E is a complex matrix while R^, are hermitian matrices. 

The NJL model is often used as a model simulating QCD. If we consider a 
single gluon exchange, it may effectively be expressed by the following four-fermion 
interaction: ' 4 ' 



^-analogue = "^(^^(^V) 



(2.4) 



where A is a (suitable) characteristic energy scale of QCD and g is the color gauge 
coupling constant. If we perform the Fierz transformation and keep only the leading 
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terms in 1/N C , this four-fermi interaction (2.4) can be rewritten into the same form 
as that in Eq.(2.1) with identification 

2 2 2 2 Mr\ 

9s=9v = 9vo = 9 • (2.5) 

( If we also keep non- leading terms in 1/N C , only the flavor-singlet vector four- 
fermi coupling is replaced by gy = (1 — 2n^/N c ) g 2 . ) We refer to this four-fermi 
interaction system with coupling relations (2.5) as 'QCD-analogue NJL model'. 

2.2. Quantization 

Let us now consider the quantization of this higher derivative system. The 
procedure has long been known since Pais and Uhlenbeck, and we here follow 
the procedure by Nakanishi. Consider generally a system which contains higher 
derivatives only in the kinetic term as follows: 

£ = ^/(^ + Ant(^), (2.6) 
where f(x) is a polynomial of the form 

n 

f(x) = aY[(x - rrij) (m 3 ;^ m k for j ^ k) (2.7) 

3=0 

and the interaction part C- m t is assumed to contain no derivatives. The well-known 
partial fraction formula 1/ f(x) = ^2j[f'(i r nj)(x — m j))~ 1 leads to an identity 

n 1 
y £ i an j [l[(x-m k )]=l = ___). (2.8) 

Using this, we can decompose the fermion field ip as 

n 

^ = ^2^3 where ipj = arjj [ (ip - m k )] ip , (2.9) 

3=0 k^j 
and then, by noting f(i$)ip = rr x {i<f) — rn^vpj independently of j, the higher 
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derivative Lagrangian (2.6) is seen to be rewritten into 



n 
j=0 

It is now clear in this form that the original higher derivative system is equivalent 
to a system consisting of ordinary (positive or negative metric*) fermion fields ipj 
with mass mj, to which the usual canonical quantization procedure is applicable. 
Then clearly the free propagator of the j-th fermion is given by ir)j/($ — mj) and 
therefore that of the original field ip is found to be given by 



F.T. <T^) = ^F.T.<T^> 



n 

I 



(2.11) 



5^-^ alU=o(y-"»i) M) 



Namely, the naive expectation that the propagator of ip is given by the inverse of 
the kinetic term is correct in this case. Conversely, if we take it for granted that 
the propagator of ip is given by i/f($), we can start from it and decompose the 
■?/>-propagator into partial fractions X]j=o ^Vj/itf ~ m j) i n every Feynman diagram. 
Then it is easy to see that the theory is equivalent to the above system (2.10) 
consisting of (n + l)-fermions ipj with mass mj. 

We can apply this general procedure to our NJL system (2.2) in various ways 
depending on which part we regard as the free kinetic term £f ree = ipf(ift)ip. We 
now discuss two ways, separately. 

A picture 

The simplest way, which we call 'A picture', is to take the original kinetic term 
+ pp/A 2 ^ as £f ree - Then f(x) = x(l — x/A)(l + x/A) and it has three 



* If we arrange the roots rrij in order of their values, the weights rjj = l/f'(m,j) take 
alternating signs and hence ipj become of positive and negative metric alternatingly. 
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roots rrij = 0, A, — A with r)j = 1,-1/2,-1/2. Accordingly, the original fermion 
field ip is decomposed into ip = ipo + + V'-A with V'o ±A denoting the component 
fermions with masses mj = 0, A, —A, respectively, and the free part lagrangian is 
written in the form: 

41 s ^(i + (212) 

= ^oip^o - 2$ A (ip - A)^a - 2^_ A ($ + A)^-a • 

This A picture hence corresponds to the following decomposition of the fermion 
propagator: 

1 11111. (2 . 13) 



The expression (2.9) of the component fermions V'O.iA m terms of the original 
fermion ip now explicitly reads 

= A _2 (^ + A 2 )^ 

1 (2.14) 
ViA^^-^i^A)^- 



B picture 

In some cases, it is more convenient to choose another form for the free kinetic 
term £f ree - Indeed when develops a nonvanishing vacuum expectation value 
(VEV) (0| S(x) |0) = ml nf 7^ 0, the fermion acquires a mass term —mipip. In such 
a case, we can take the following lagrangian as the free kinetic term by including 
the mass term: 

Then, for this choice, we have 

f(x) = x(l - x 2 /A 2 ) -m = -A~ 2 (x - m Q )(x - m x ){x - m 2 ) , (2.16) 
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and the kinetic term (2.15) is rewritten into 



£free = % Vo(«# - ™o)^o + Vi - m^ip! + r] 2 1 ip 2 (ift - m 2 )4>2 ■ (2.17) 

We refer to this as 'B picture', which corresponds to the following decomposition 
of the fermion propagator: 

1 111 . 

T, ^T = ^° J +r li J + r ?2 j ■ 2 - 18 

m-#f(l-^) rn -p mi-f m 2 - f 

An inconvenience for analytic treatment in this choice is that we have no simple 
explicit expressions for the three masses rrij (and hence the weights r)j = [f'irrij)]^ 1 
also); they are determined by Eq.(2.16) and their explicit expressions can only 
be given by a complicated Cardano's formula. However, for the purpose to do 
practical calculations, we can make use of computer. Then such complication is of 
no problem and the decomposition (2.17) reducing the higher derivative system into 
the usual fermion system provides us with very efficient tool to calculate various 
quantities, as we shall see explicitly in Sect. 3. 

We assume henceforth that we are labeling the three masses (mo, mi, 7722) in 
order such that they approach to (0, +A, — A), respectively, as the mass parameter 
m goes to zero. So ipo with mass mo is the physical fermion with positive metric 
and the other ipi ;2 are unphysical fermions with negative metric. We should note 
the fact that possible value of the physical fermion mass mo is bounded from above 
by 

m < = 0.57735 • • • A, (2.19) 
v3 

although the mass parameter m can be arbitrarily large in principle. This is because 
the root mo becomes complex beyond this limit. Indeed, the cubic polynomial f(x) 
in (2.16) has two stationary points at ±A/\/3 and the root mo, if being real, has 
to lie in between them. This bound (2.19) is not a defect of the present definition 
of the NJL model but is of physical significance. The NJL model is defined with 
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an ultraviolet cutoff in any case and can only describe physics below it. Therefore 
the physical fermion mass mo generated by the NJL model dynamics can be at 
most of the same order as the cutoff and otherwise becomes unreliable. 

The expression (2.9) of the component fermions ipj in terms of the original one 
ip now reads 



f^o\ 




/ Vo 


rjomo 


in 








w 




\V2 


m^2 



(2.20) 



V V / 



where we have used the following relations for the three masses rrij which directly 
follow from the defining equation (2.16): 

mo + mi + rri2 = 0, mQirii + m\rri2 + rri2m§ = —A 2 , m§m\rri2 = —mA 2 . 

(2.21) 

For later convenience, we here cite some formulas for the masses rrij and the weights 
r]j. Using the relations (2.21) and m 3 = A 2 {rrij — m) for Vj, one can easily derive 



3=0 

£ 

%>j 



2 2 



A 4 









X>5 = 


-3mA 2 , 


X>4 = 2A' 


3=0 




3=0 


m 
rrij 


9 
3 


etc.. 



(2.22) 



The defining equation of the weights rjj together with these mass relations lead to 
the following relations: 



/^n_ 2 m 1 
\rrij A 2 



2 2 

3=0 3=0 
2 2 2^2^ 

V ^m 2 = -A 2 , V ^m 3 = 0, V -^L = -, J2^2=— , etc.. 

j=0 j=0 j=0 3 3=0 J 

(2.23) 
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2.3. Current Operators 

We have emphasized that our definition of NJL model respects the chiral sym- 
metry. It is indeed clear that our Lagrangian (2.2) is invariant under the chiral 
U(nf)ji x U{ri{)i J transformation. Let us derive current operators corresponding to 
this symmetry since they take different forms from the usual ones because of the 
presence of higher derivatives. 

To do this, the simplest way is to use the A picture. Then the Lagrangian is 
given by (2.12) (plus interaction term containing no derivatives) to which the usual 
Noether procedure applies. In view of the expression (2.14) for the component 
fields V'OiA m terms of the original ip, we see that the vector transformation 
tfj — > ip' = e ip (6 = Yla @ a ^ a ) on the original fermion field ip is realized on 
the component fields VfyiA m the same form, ^o,±A - ► V'o ±A = e yO,±A, while 
the chiral transformation tjj — > if}' — e^^ip is realized on the component fields a 
bit differently as follows: 

A + ^-A = e i75e (^A + ^-A) (2.24) 
A-^-A = e-^(^A-^-A) • 

Having found the transformation law of the component fields ^o.iAj we can apply 
the usual Noether procedure and obtain the corresponding currents. The vector 
current is found to be 

f a = t/SoAVVo - 2^aA°7^A - 2^_ A A a 7^-A , (2.25) 
and the axial current is given by 

& = ^ A a 7 M 75^o - 2^ A A a 7^75^-A - 2^_ A A Vts^A • (2.26) 
We can rewrite these current expressions in terms of the original fermion field ip 
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by using relations (2.14) and find 



JL = <M V75^ + ^ W • AV75^ + <A P A V75^ + Vt5D V-a] • 

(2.27) 

Although not being well-known, there is in fact a direct and general proce- 
dure for deriving the Noether current (without decomposing ip) for generic higher 
derivative systems, which we show in the Appendix. We can see that the Noether 
current obtained by the direct procedure coincides with the above (2.27) up to 
so-called ambiguity term of the form d u f^ u with an antisymmetric tensor f^ v . 

In some cases we need current expression written in terms of the component 
fields ipj in B picture. This can be obtained if we invert the relation (2.20) between 
ipj and (ip, jzV, □?/>): 







/ml 


9 

m\ 


m\\ 








= rriQ 


m\ 


rri2 




V 


i> / 


\ 1 


1 


1 / 





(2.28) 



Substituting this into (2.27), we find the currents in the B picture: 



j=0 

2 2 



2m,'m, 



3=0 



i,j=0 



(2.29) 
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2.4. Coupling to External Gauge Fields 

It is easy to couple external gauge fields to the system by gauging the chiral 
U(nf)-R x C/(rif)L symmetry. This is achieved simply by replacing the kinetic term 
by the covariant derivative one: 

(2.30) 

where = V^A a and A^ = ^ a *<4^A a are the vector and axial- vector gauge 
fields of the chiral U(rif)nX U(nf)i group. If we gauge a part of the group, the gauge 
fields should also be understood so; for instance, if we couple only the photon to 
the system, then V M = eQA^ and A^ = with Q being the charge quantum number 
matrix of the fermion. Note that the covariant derivative of course depends on 
the operand and D^ip = (d^-iV^-iA^^ip while D^Ipijj = (d^-iV^+iA^^Ipip. 

Here two remarks are in order. First we note that the above covariant kinetic 
term contains the terms linear in the gauge fields in the form, 

+ A "A + ( total derivatives), 

and the currents ja and j§ just coincide with the above ones given in (2.27). 
Second, in the above covariantization we started with the expression ipift(l + 
(<f)ft/ A 2 ))ip. But if we started with an equivalent one + (□/A 2 ))-?/', we 

would have obtained ^%Jp{\ + (D 2 / A 2 ))?/>, which is different from the above one 
^>iJp{\ + {JpJp I A 2 ))-?/>. The difference is, however, seen to be a non-minimal inter- 
action term (Pauli term) of the form like F^^^^f^) with a field strength F^ V) and 
hence the currents defined by the linear terms in the gauge fields again coincides 
with the above ones up to ambiguity terms of the form d v f^ v . 

We have seen that our definition of the NJL model by higher derivative is 
made consistent with the gauge symmetry also. However, we should note the 
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fact that the vertex functions of the external gauge fields alone are not yet well- 
regularized by our higher derivative fermion kinetic term. Consider generally a 
fermion one-loop diagram which contains n vertices of 'mesons' E, and and 
m vertices of external gauge fields V M and A^. The fermion propagator behaves 
as ~ A;~ 3 at high loop momentum k and the vertex factor for 'mesons' contains no 
momentum. But the point is that the vertex of the external gauge fields contains 
second powers of fermion momenta as the currents (2.27) show. Therefore the 
diagram has superficial degree of divergence 

uj = 4 — 3(n + m) + 2m = 4 — 3n — m. 

This is non-negative when (n — l,m = 0,1) and (n — 0,m = 1,2,3,4) aside 
from an irrelevant vacuum graph case (n — m — 0). The linear or logarithmic 
divergences for the former cases (n — l,m = 0, 1) are in fact absent because of 
Lorentz, chiral and gauge invariances. Indeed, first, the (n — l,m — 0) case cor- 
responds to 'meson' tadpoles (1-point functions), and E-tadpole vanishes by the 
chiral symmetry while vector and tadpoles vanish by Lorentz invariance. 
Second, the (n — l,m — 1) case corresponds to 'meson'-gauge transition 2-point 
vertex functions; when the meson is E, the 2-point vertex E-V M or T.-A^ van- 
ishes by chiral symmetry, and when the meson is R v or L v , the 2-point vertex 
like (R + L) v -V^ should be proportional to (g^p 2 — VyPv) by gauge invariance. 
Therefore they, being superficially logarithmically divergent, become convergent 
actually. The divergences thus occur only when (n — 0, m — 1, 2, 3, 4), i.e., for the 
diagrams consisting of external gauge fields alone. For those diagrams the present 
higher derivative kinetic term does not improve the divergence situation at all; 
however higher the fermion kinetic term derivatives is chosen, the gauge boson 

vertex also gets to contain higher derivatives accordingly. This is a well-known 

[7] 

fact in the higher derivative regularization. We can however regularize those di- 
vergence using dimensional or Pauli-Villars or any other suitable regularization* 

* Note that this is a mere regularization which is to be removed eventually by the usual 
rcnormalization. 
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For definiteness, we adopt dimensional regularization in this paper. 



3. Example Calculations 



3.1. Effective Action and Potential 

As we have seen in the above, the naive Feynman rule is correct for this type of 
higher derivative system. This also implies that the usual Feynman path integral 
expression for the Green function generating functional is valid for our system with 
lagrangian C given by (2.2). Therefore the effective action T for the 'meson' fields 
S, £t, R^ and is given in the leading order in 1/N C simply by integrating over 
the fermion field if) (and if)): 



derstood that this vector singlet term is always accompanying the vector term 
(A 2 /2) tr( J R^ J R /i + L^L^)] Precisely speaking, the second TrLn term standing for 
the fermion one-loop diagrams is not made convergent yet by our third order deriva- 
tive propagator. But if we expand it with respect to A4, the ultraviolet divergence 
appears only in the first two terms; the zero-th and first order terms in M.. Since 
the first order term vanishes by chiral symmetry and Lorentz invariance as we ex- 
plained above, the divergence in fact appears only in the zero-th order term which 
is a field independent constant we can throw away freely. 

The effective potential V(E), whose stationary point determines the VEV of 




(3.1) 



+ -^TrLn[^(l + nL)-M]. 
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£, is easily found from (3.1) to be given by 

^ tr(EtE) _ 2Arc |_g ltrln[E , E _ p2(l _^ ] , 



(3.2) 



where tr in the second integral expression denotes the trace only over the flavor 
space. The scalar field £ is generally expanded into 



a=0 a=0 

and we see that the VEV of E is always brought into a real and diagonal matrix 
form by the chiral £/(rif)R x f(^f)L rotation and further that the diagonal matrix 
is in fact proportional to unit matrix since each of the diagonal values is separately 
determined by minimizing exactly the same form of potential. Therefore, without 
of loss of generality, we can substitute 

(E) = a°X° = ml^ (3.4) 

into the effective potential, where m is a mass parameter which will yield a fermion 
mass term —rmpi/j. Then we find the effective potential to be 



1„ (E = ml „,) = |L m2 _ 2Nc J Jg. ln [m » _ p2(l _ . (3.5) 



The stationary condition of the potential leads to the following gap equation which 
determines the non-zero mass value m corresponding to spontaneous chiral sym- 
metry breaking: 

A 2 f d 4 n 1 

1 dP ' 0. (3.6) 



4glN c J ? (2vr) 4 m 2_ p 2( 1 _g) 2 



Up to here the story is formally the same as in the usual NJL model case and 
the difference is only in the higher derivative term in the fermion propagator. As 
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noted in the previous section, the higher derivative propagator can be reduced to 
the sum of ordinary fermion propagators and then we can evaluate the quantities 
in the same manner as usual. For instance, consider the effective potential in (3.5). 
If we use Eq.(2.16) with x = ±$ substituted, we can evaluate it as 

I A 2 ^ 

— V(E = mln) = -^m 2 + N C J2 *bK) , (3.7) 

where rrij are the three masses of the component fermion ipj, determined by 
Eq.(2.16), and the function Fo(mj) is the usual fermion vacuum energy given by 



— — — tr \n(mj — 4) = — 2 / — — — \n(m 2 - — p 2 ) 
i(2vr) n K J rj J i(2ir) n v 3 F ' 



I , /, . , 3 



(3.8) 



4 

r m 



L + ln-|-- + 0(e), 



16tt 2 3 \ A 2 
with L standing for the 'divergent' part: 

9 1 

L = In A — -, 

e 

1 1 , , x 4-n 

- = 7 + ln(47r), e = — - — , 7 : Euler constant . 

6 6 Zi 



(3.9) 



Here note that we have used dimensional regularization to evaluate the contribu- 
tions from the component fermions. This is always necessary since those contribu- 
tions are separately ultraviolet-divergent, although the sum is generally convergent. 
In this case of effective potential, however, the sum still contains an m-independent 
divergence as mentioned above. In fact, the sum Y2,j=Q^\ m j) contains the diver- 
gent part — (L/167T 2 ) X^ m j an d it is seen to be m-independent constant because 
of an equality Yl j m j = 2A 4 in (2.22). Discarding the m-independent divergence 
and constants, we find the potential (3.7) to yield 



— V(E = ml nt ) = -^m 2 - Vmjln^- . (3.10) 
n f v l ' 2ql 16tt 2 ^ 3 A 2 



The constant is adjusted so that this becomes zero as m — > in which (mo, mi, m.2) 
— > (0, +A, —A). The gap equation (3.6) can also be rewritten as follows if we use 
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an equality [m 2 — p 2 (l — (p 2 /A 2 )) ] = J2jVj( m j/ m )[ m 'j ~ P 2 ]" 1 which follows 
from taking trace of both sides of Eq.(2.18): 

A 2 1 2 m 3 - m 2 

—5 — = o> Vi— n-4 . 3.11 

Ag 2 N c 16vr 2 ^ lj m A 2 V ; 



3.2. Two-point Functions 

The same technique applies also to other quantities. For instance, the n-point 
Green functions on the spontaneously broken vacuum is obtained by expanding 
the effective action (3.1) around M. = m. So the fermion propagator is given by 
i[$f(l —p 2 / A 2 ) — m] _1 which can be decomposed into the same component fermion 
propagators as above. Let us demonstrate this by calculating the 2-point functions 

(2) 

To- or 7r of scalar a and pseudoscalar it on this vacuum. For every flavor components 
a a and ir a , independently of a = 0, 1, • • • , n 2 — 1, they are given by 



r ( 2 ) ( v 2 )--^ 
Kr 2g 2 

_Nc 
2 



d 4 k 
i(2ir)< 



tr 



1 

il5 



m 



A 



1 

il5 



m 



A 



(3.12) 



where $]a is a shorthand notation for #(l — k 2 /A. 2 ). The second term standing 
for the fermion one-loop diagram can be evaluated by using the B picture decom- 
position (2.18) for each propagator as follows: 



second term 



N c JU f d n k 



i,j=0 
2 



1 



i,j=0 



«75 J m - k* { % J mj - + |f) . 

G 

2. 



^2 2^ ^ [K + m i ± m * m ^ - y) L " ( — 2 -) 



+ F±(p ]m h mj) 



(3.13) 
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with 



/' 
dx (2mfj ± mirrij — 3x(l — x)p 2 ) In • 



mjj —x(l — x)p 2 



A 2 



\ — x)m 2 + xm 2 . 



(3.14) 

All the divergent terms proportional to L in (3.13) are seen to vanish if we use the 
?j=oVj = and J2 2 j= 



identities Yl^oVj = an d Y^j=Q r lj m 3 = i n (2.23) and therefore we obtain 



r (2) (v 2 ) 



A 2 N, 



+ FT ViVjF±(p 2 ;m,mj) . 



(3.15) 



2^ 8tt^ .f 

yb 1,3=0 

In quite the same way we can calculate 2-point functions for the vector and axial 
vector mesons, = (R^ + and = (R^ - assuming g% = g^ Q , 

they are given independently of the flavor by 

r(2) ^(p 2 ) 



{5} 



9 g\ 2 2-, ^ 7 *(2tt)4 [ \ 7 m 75 J m , _ y \ ri5 j mj - [# + f) \ 



9 



(lis 



A 2 _ 
9 2 v 8vr 2 



A", 



(3.16) 



where 



Gj^j (p ; mi, rrij) = j dx (m^ =f mirrij) In 



mfj — x(l — x)p 2 



A 2 



iJ(p ;mi,rrij) = / cix 2x(l — x) In 



m?- - x(l - rc)p 2 
A 2 



(3.17) 



Although these Eqs.(3.15) and (3.16) are not so explicit expressions since the roots 
rrij are implicit functions of m and A, they give closed expressions for the 2-point 
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functions which allows straightforward numerical evaluations. (In particular, the 
parameter integrals over x in (3.14) and (3.17) can be carried out analytically and 
the functions F±, G ry-i and H are given explicitly by elementary functions.) 

w 

We here note some points on these 2-point functions. Firstly the pion is 
indeed massless, Y^\p 2 = 0) = 0, as it should be because of chiral symme- 
try. This can immediately be seen by putting p = directly in (3.12) and 
comparing it with the gap equation (3.6). [It is more complicated to see this 
if we use the expression (3.15), although being possible by the help of a for- 
mula J2jVj/( m j + m i) = l/2m (mi-independent) following from the relation 
J2jVj/( m j — x) — — A 2 /nj( m i — x )-] Thus the pion 2-point function behaves 

as vf\p 2 ) = Z~ 1 p 2 + 0(p 4 ) around p 2 = and the coefficient Z~ x defines the 

i 

(inverse of) wave-function renormalization factor of our n field as 7T = Z£ir r . Since 
there is no genuine Yukawa vertex correction in the 1/N C leading order, this wave- 
function renormalization factor also determines the Yukawa coupling g^J,^ defined 

_ i 
by Ant = ~9Tnpip^ / y5 7T r^', namely, g^^ = Z% . Evaluating the coefficient of p 2 in 

(3.15), we find 



7-1 



i,j=0 



(3.18) 



dp 2 



(0;rrii,mj) 



-mf(mi + 2m 



2{mi + m,-) 3 (mj — m,- 



+ (i <-> j) 



+ 



mirrin 



+ 



2(m, + rrij) 2 12 



2 A 2 6 



(3.19) 



Secondly, an interesting point in our definition of NJL model is that the scalar 
a meson develops a pole below the two fermion threshold 2mo- Namely, the famous 
NJL mass relation m G = 2mo no longer holds here. We can find the a meson mass 

(2) 

numerically as a zero of To- in (3.15) and show the result in Fig.l. Note that 
the NJL mass relation holds only in the limit mo/ A — > 0. [The approach to the 
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Fig.l. Meson mass squares ml vs fermion mass mo- er, V and 
A denote scalar, vector and axial-vector meson boundstates, 
respectively, (mo /A is bounded by l/y/3.) 



Finally, as for the vector and axial- vector meson 2-point functions (3.16), we 
note that the one-loop contribution parts do not have gauge-invariant form oc 
(d^P 2 ~V[iVv) even in the limit m — > contrary to the true gauge fields. This is of 
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course because the apparent 'gauge-invariance' for those mesons is violated in the 
higher derivative terms in the present model. We have also plotted in Fig.l the 
masses of the vector and axial-vector mesons determined as zeros of in (3.16) 
for the case (?| = corresponding to QCD analogue model. It may be of some 
interest to note that, if we take A = lGeV, 900MeV, 800MeV and f„ = 93MeV as 
inputs, the present NJL model gives the following values for the fermion mass mo 
and the fermion pair VEV: 

A = 1 GeV : m = 250 MeV, ( - <^> 1GeV ) 1/3 = 253 MeV 
A = 900 MeV : m = 273 MeV, ( - (V^> 1GeV ) 1/3 = 245 MeV (3.20) 
A = 800 MeV: m = 315MeV, ( - <^> 1GeV ) V3 = 237MeV , 

where Eq.(4.7) below for f n and Eqs.(3.18), (3.6) and (2.3) are used. Eq.(2.3) 
gives — (V^) — (^■ 2 /fi's) m ' which we regard as the fermion pair VEV renormal- 
ized at the cutoff scale A. The fermion pair VEV's cited here are renormalized 
ones at /i — 1 GeV which we calculated using the following formula although the 
renormalization effects are small: 

(bKAVA^,,) J 7 s m (3 ' 21) 

with Aqcd — 500 MeV, where 4/9 is the anomalous dimension of tpvp in three 
flavored 577(3) QCD. The empirical value for ( — (^) 1Ge v) ^° he compared 
with these is 225 ± 25 MeV. 
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Fig. 2. Feynman diagram contributing to /, 
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Note that the RHS of Eq.(4.2) should be evaluated at the mass shell p 2 = 
in order to give the decay constant f^. But the RHS itself, corresponding to the 
diagram in Fig.2, is defined even off the mass shell and hence the Eq.(4.2) may be 
regarded as defining an 'off-shell' decay constant f n (p 2 ) which coincides with the 
true f n at p 2 = 0. With this understanding, let p^ be off the mass shell for a while 
and multiply both sides of Eq.(4.2) by ip^. Then we can obtain 



Up 2 )p 2 = ^ 



_<£k_ 

i(2ir) 1 



-tr 



m — 



+ 2m tr 



1 

]a m-[ 
1 



A 



1 



-75 



m-[#]A m-[y + j]A 



75 



(4.4) 



(4.5) 



by using an algebraic identity 

p^A^h, k + p) = $ + ^]a75 + 75$]a 

= (m - [#]a)75 + 75 ("i - W +^]a) - 2m75 . 

Note that this identity (4.5) is just a higher derivative case version of the usual 
tree level Ward-Takahashi identity for the axial vector current. The two terms 
in the first trace in (4.4) become equal to each other by shifting the loop mo- 
mentum (which is allowed now) and they give — 2m(A 2 / g^N c ) x (g^^N c /2) = 
—mg^^k 2 jg\ owing to the gap equation (3.6). So we find 



fn{p 2 )p 2 = 2m g n 



A2 

2d 



i(2it) 1 



tr 



m 



75" 



m 



75 



(4.6) 

But the quantity in the curly bracket here is just the same as the 2-point function 

T^\p 2 ) in Eq. (3.12) of pseudoscalar ir, which we know behaves as Z~ l p 2 + 0(p 4 ) 

i 

around p 2 = 0. Therefore, taking also account of the relation g^^ = Z% , we find 

2m 



(4.7) 



Although we have derived this relation (4.7) by a rather explicit computa- 
tion in the leading order in 1/N C , it is in fact a direct consequence of the chiral 
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symmetry alone. Indeed, using canonical commutation in A (or B) picture and ex- 
pressing n and a fields in terms of fermion field, we can derive the chiral symmetry 
transformation law of the NG boson field n b : 

6(x )[ij 5a (x), rr b (0)] = 6\x)tr [a & (AV(0) + a(0)A a )] . (4.8) 

Then, by using VEV of this equation and (0| a(x) |0) = ml nf , we see that the 
current conservation d^j^x) = leads to 

^ <0| T j^(x)n b (0) |0) = 6(x°) <0| [ij° 5a (x), vr & (0)] |0) = Am8\x)8 b a . (4.9) 

This implies that the 2-point function (0| T j§ a 7r b |0) in momentum space contains 
a massless pole and is given by (p^ / 'p 2 )Am5 b a . But, on the other hand, such massless 
pole comes from the NG boson intermediate state and hence the pole residue should 
be 2/ 7r ]9 /i ■ Z£8 b . We thus get an equality 2f n Z^ 2 = 4m, the same relation as the 
Eq.(4.7). 

4.2. 7T° — > 27 Amplitude 

We now turn to calculate the amplitude of 7r° — > 27 decay. The neutral pion 7r° 
1 

corresponds to Z£% s in the present notation and the photon 7 couples to (e times) 
the vector current ja(x) with the flavor matrix A a replaced by the quark charge 
matrix Q = diag(2/3, — 1/3, • • •). The simplest way for calculating this amplitude 
is to use the B picture since the vector current is diagonal with respect to the index 
j of the component fermions ipj even in the B picture. The amplitude is calculated 
by a triangle diagram and is given by 

A 3 2 

= ~e 2 g^N c tr (QQ-) ^ VjT^ip, q; mj) , 

j=Q 



T^(p,q; mj ) 



r d n k 






[■ 



7" ~~Ti^ 7u 1 j\ (~ j 75) 



m j m j -H m j-(ty + f) 



(4.10) 
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The evaluation of the integral of the amplitude T^(p,q;m,j) is well known;' 10 ' it 



gives 



T^(p, q- mj) = -^—^q^p , (4.H) 



so that we find 



3=0 J (4.12) 



Here we have used J2j r ]j/ m j — l/ m i n Eq.(2.23). If g ir ^/2m = 1/ f n , this exactly 
reproduces the well-known low energy theorem for 7r° — > 2-y. This is indeed the 
case because of the relation (4.7) as we confirmed in the above. 

5. Simplified Treatment for the Case A ^> m 

[91 

In some applications of NJL model (as in the top condensation scenario, for 
instance), the cutoff A is much larger than the scale of the mass m and momenta 
which we discuss. In such cases we need not keep terms which are suppressed by 
1/A 2 in the effective action. Then we can take A 2 — > oo limit in all the terms not 
diverging in that limit. This simplifies the calculations considerably and the cal- 
culated results become very explicit ones containing no longer the implicit masses 
mj. 

We now show how to evaluate the effective action (3.1) in such cases. Only 
problem is the fermion one-loop term, which we expand as follows using the nota- 
tion [ip] A = ip(l + pp/A 2 ): 



TrLn[[ij% - M] = Tr 



Ln[[ifl A ]--^---( M 




(5.1) 
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The first term in the RHS is an irrelevant constant and the second term vanishes as 
explained before. Now we consider the A — > oo limit of this quantity (5.1). Recall 
the propagator decomposition in A picture: 



and call the LHS regularized propagator, the first term in the RHS unreg- 
ularized propagator and the second and third terms with masses ±A regulator 
propagators. In the last term Tr[ — ^ n>5 (l/n)(7W/i[^]A) n ] in Eq.(5.1), we substi- 
tute (5.2) for each propagator factor l/i[ip]\ and then it becomes a sum of various 
terms consisting of the unregularized and regulator propagators. But, from dimen- 
sion counting, each term is ultraviolet convergent and therefore we can take the 
A — > oo limit directly for the loop integrands (i.e., inside the functional trace Tr). 
Clearly then all the terms containing the regulator propagators at least once drop 
out and only the term Tr[ — ^2 n >5(^/n)(A4/ip) n ] survives. If we add to the Tr 
operand of this term the quantity 



1 111 11 



(5.2) 



[ip] A ip 2ip-A 2ip + A ' 




then it reproduces the one-loop term TrLn \ip — JA\ of the unregularized fermion. 
We thus find that 




up to an irrelevant constant. [Precisely speaking, the last term and the second last 
term Tr[ — (l/A)(M/i[p]A) 4 ] in Eq.(5.1) separately have an infrared divergence 
and the discussion above is not so rigorous. But actually a more careful argument 



can justify this expression (5.3) as is inferred from the fact that in this expression 
(5.3) the infrared divergence cancels between the terms Tr[— (l/4)(A1/i[^]A) 4 ] 
and Tr[(l/A)(M/ip) 4 ].} 



The Tr operation (loop integral) of this expression (5.3) as a whole of course 
gives a convergent quantity. But the integral for each term separately is divergent. 
In practice, however, we can evaluate the integral for each term separately if we 
use dimensional regularization. The convergence as a whole implies that the pole 
terms 1/e appearing from each term should cancel eventually. We can use this fact 
to check the calculations. 

We can evaluate the 'counterterm' r count defined in Eq.(5.3) in a closed form 
if we discard 0(1/A 2 ) terms since then only dimension 2 and 4 operators survive. 
By using Eq.(5.3), the effective action (3.1) reads: 



Performing a straightforward (but a bit laborious) calculation we find that the 
counterterm r coU nt is given by 




(5.4) 



+ -^TrLn[^ — M] + T count + 0(-^) . 
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count — -L count ~<~ 1 coun t 1 count ' 



.(2) 
count 



16tt 2 



2A 2 tr(Sts) + (L - ^) tr(^Et^E) - ^A 2 tr(#% + 



6 tr 



(^i? u ) 2 l + (R - L) 



.(3) 
count 

Nr. 



16lY 2 



,(4) 
count 



16tt 2 



i(L - -) tr [i^(£t + L^(Sa M St) 



+ -(L - 2) tr (c^iT - cW) i2„] + (i2 — > L) 



_ (L - ^) tr [(S+S) 2 ] + (L — ^) tr (i? 2 £t£ + L 2 £Et) 
17 



2(L--)tr 



' ' (5-5) 

where T^ nt with n = 2,3,4 stand for the contributions from (l/ra)TrLn[ — 
(M./[i<f)}A) n + (M/ip) n ], respectively, and L is the divergence factor In A 2 — 1/e 
introduced in (3.9). It should be noted that the divergent terms proportional to L 
in r coun t are combined to yield just the following 'gauge covariant' form: 



1 



N c - 
16tt 2 " 



tr(^St^S) --triF^F* + F L ^F^) 



(5.6) 



with 



Z^E = c^E - iL^E + *Ei? u 



F* = d^X v - d v X^ - i[Xp, X v \ {X = R,L). 



(5.7) 



This 'gauge covariant' form for the divergent parts is a reflection of the fact that the 
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fermion one-loop term oc TrLn[ — Ai] formally reduces to a gauge invariant 
form TrLn[r^(<9 M - iR^Vr - iL^V L ) - T,V R - T^V h ] in the A -> oo limit as if 
and are the gauge fields of the local U (rif )r x U (rif )l symmetry. The finite 
parts, however, do not have such gauge invariance, of course. 

This form of the effective action (5.4) with counterterm r coU nt is a very conve- 
nient one. The effect of the cutoff by higher derivative term is now isolated in the 
counterterm r count and it is given explicitly in (5.5). All we have to do then is to 
calculate the simple fermion one-loop action TrLn[ij$ — Ai ] with no cutoff by using 
dimensional regularization. The divergences appearing there will be automatically 
canceled by the contributions from the counterterm r coU nt- 

Let us demonstrate the simplicity of this way of calculations by taking an 
example — the effective potential V(m) = l / (S = ml nf ). The effective potential 
contribution from the unregularized fermion one-loop is easily evaluated as 

1 f d n p 

— Kmrcg.l-loopM = -Ac / trln[m - f[ 

f d n p , r n 9l A^ ,/ 3 m 2 \ 
= - 2N <Ji(2^ hl [m " P 1 = "I6^ m ( L - 2 + A*) 



(5.S 



while the contribution from the counterterm is found by substituting E = ml ni 
into (5.5) to be 

^count(-) = -^(2- 2 A 2 - - 4 (L " y)) • (5.9) 

We see the divergence oc L in (5.8) is actually canceled by this counterterm. Adding 
the tree contribution nf(m 2 A 2 /2#§) also, we have 

1 , s A 2 9 A" c / 9l9 A , m 2 4 N \ 

—V(m) = -^m 2 % ( 2m A + m 4 ( In — T + -) J . 5.10 

n f V ; 2q 2 16n 2 \ V A 2 3V V ; 



The validity of this expression is of course limited in the region m /A 2 < 1. We 
draw in Fig. 3 pictures of this potential and the 'exact' potential (3.10) in Sect. 3, 
for comparison. 
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V(m) 




Fig. 3. Effective potentials for the case gs = 1.0193(?g r above the 
critical coupling g§ r = 2-7T / 'y/N^. Vertical and horizontal axes 
give scaled ones: V(m) x (16tt 2 • 10 3 /iV c A 4 ) and m x (10/A). 
Broken line denotes the approximate one (5.10), solid line the 
'exact' one (3.10) and dotted-broken line the one (5.11) of the 
cutoff theory. 

It may also be of some interest to compare this potential (5.10) with the usual 
effective potential obtained by the simple cutoff of the loop momentum which reads 




We have also drawn this potential function in Fig. 3. [Its deviation from the present 
'exact' one (3.10) in the large m region in Fig.3 reflects the difference in the meaning 
of the cutoff A in both theories.] 
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6. Conclusion 



In this paper we have defined the NJL model by introducing higher derivative 
fermion kinetic term. We have clarified some basic aspects of the theory such as 
quantization, current operators etc., and developed two calculation methods which 
make the evaluation of the diagrams in this higher derivative system no more 
difficult than in the usual first order derivative case. 

We have emphasized that the present formulation of the NJL model suffers 
from no ambiguities in the loop momentum assignments and keeps the important 
chiral and gauge symmetries. We explicitly demonstrated this by calculating the 
7iq — > 27 amplitude and confirming that the low energy theorem holds. 

We restricted ourselves in this paper to the system with exact chiral symmetry. 
But there will be no problem in extending the present formalism to the cases where 
the fermions have explicit chiral symmetry breaking masses. When a fermion has 
an explicit mass M, we think it best to take the kinetic term in the form: 



Indeed, then the fermion really has mass M in the absence of interactions and 
satisfies correct normalization. The fermion momentum is effectively cut off as 
\p 2 — M 2 1 ^ A 2 around the mass shell. This form will be suitable also for describing 
heavy quarks which are much studied recently in connection with heavy quark 
symmetry! 4 ^ 




(6.1) 
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APPENDIX Noether Current for 
a generic Higher Derivative System 

We consider a generic system whose action contains arbitrary order derivatives: 

S(<j>) = J d 4 x£((j), <9 M 0, d^(f>, d^ p (p, •••) , (A.l) 

where <fi stands for a collection of fields (whose index is suppressed) and we use 
abbreviations like 



(A.2) 

weight 1 



The suffix 'weight 1' in the latter means that we keep always the weight to be 
one irrespectively of whether the n indices fi\,fi2, • • • , take the same values or 
not; namely, for the case C = a^d^cj), for instance, dC/d{d\\<p) = a 11 and 
dC/d(d 12 <p) = a 12 + a 21 , but «9£/<9(<V0)| weightl = (aT + a"")/2! always. The 
functional derivative of the action S with respect to is given by 



^ = ^-fy(<f£^)+M*£ ; H- d^ p (5C^) + --- , (A.3) 



and the Euler-Lagrange equation of motion is written as 5S/5(p = 0. If we perform 
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an infinitesimal transformation 



i>, the Lagrangian C changes as 



8C 

SC = —50 + (SC ; ^50 + (SC ; H^50 + (SC ^ p )d^ + ■■■ 

= ^50 + [(SC ; %50 + d^SC >n ■ 50] + [(SC ; n<V^> " 9^(SC ; H ■ 50] 
+ [(SC ^ p )d^ p s<t> + d^ p (SC -m) ■ s<p] + ■ ■ ■ . 

(A.4) 

To rewrite this, we introduce a generalized 'both-side' derivative defined by 



F d fiiH2—(in G — Fd piP2 ... Pn G d pi F ■ d P2 ... Pn G 



+ 9jj, lfJ , 2 F ■ 9 M3 ... Mn G h (-Td^^.-.^F ■ G 



(A.5) 



for arbitrary two functions F and G. This derivative is no longer symmetric under 
permutation of the indices but enjoys an identity 

dii [ F ^~^d av .. an G] = F^- a "d fiai ... an G + (-) n d fiai ... an F^- a - ■ G (A.6) 

for functions F^ ai "' an totally symmetric with respect to the n + 1 indices /i, ai, 
• • • , a n . Applying this identity we can write (A.4) as 

5£ = —50 + <9 M [(5£ *)6<i>] + <9 M [(5£ ^) dj<f>] + <9 M [(5£ d„ p 50] +■■■ 



■£(6C'->"»-<>")d ai ... an 6<t> 



.71=0 



(A.7) 



If the transformation — > + 50 with 50 = e a 5 a (e a : x-independent transfor- 
mation parameters) leaves the lagrangian invariant, SC = 0, then we obtain from 
this an identity: 



^ = ^(5£^—)Sa,..a„5a0. 



(A.8) 



n=0 



This ja gives a generalized Noether current and is seen to be conserved if the 
Euler-Lagrange equation SS/ 50 = is satisfied. 
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For completeness we show here that this generalized Noether current coincides 
with the source current of the gauge field if it is introduced by the usual covari- 
antization procedure replacing the derivative by a covariant one: — > = 
—iA a ^T a for the case 5 a (j) = iT a (p with a certain representation matrix T a . Then 
noting 



D^D^ ■ ■ ■ D^ n <p = d^ lfl2 ...^ n 



k=l 



(A.9) 

we find the term linear in the gauge field in the covariantized lagrangian £ cov (0, A) 
is given by 



n=0 fc=0 



n=0 



.k=0 



+ (tot. der.) , 



(A.10) 

where we have performed a 'partial integration' in going to the second line and (tot. 
der.) denotes a total derivative term appearing then. We note that the quantity 
multiplied by A^ in the first term is just identical with the above Noether current 
and so we have shown 

5S C0V [<f>,A] 



5A* 



A=0 



J a 



(A.11) 
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